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Abstract 

Neutrino energy losses through neutral weak currents in the triplet-spin superfluid neutron liquid 
are studied for the case of condensate involving several magnetic quantum numbers. Low-energy 
excitations of the multicomponent condensate in the timelike domain of the energy and momentum 
are analyzed. Along with the well-known excitations in the form of broken Cooper pairs, the 
theoretical analysis predicts the existence of collective waves of spin density at very low energy. 
Because of a rather small excitation energy of spin waves, their decay leads to a substantial neutrino 
emission at the lowest temperatures, when all other mechanisms of neutrino energy loss are killed 
by a superfluidity. Neutrino energy losses caused by the pair recombination and spin-wave decays 
are examined in all of the multicomponent phases that might represent the ground state of the 
condensate, according to modern theories, and for the case when a phase transition occurs in the 
condensate at some temperature. Our estimate predicts a sharp increase in the neutrino energy 
losses followed by a decrease, along with a decrease in the temperature that takes place more 
rapidly than it would without the phase transition. We demonstrate the important role of the 
neutrino radiation caused by the decay of spin waves in the cooling of neutron stars. 
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I. INTRODUCTION 



Usually neutron stars consist mostly of a superdense neutron matter which is in /3 equi- 
librium with a small fraction of protons and contains the triplet-correlated superfluid con- 



densate of neutrons below some critical temperature . For a long time, it has been 

generally accepted that the pair condensation in the superdense neutron matter occurs into 
the state (with a small admixture of ^^2) with a preferred magnetic quantum number 
nij = 0. This model has been conventionally used for estimates of neutrino energy losses in 
the minimal cooling scenarios of neutron stars 

a, a 

. During the last decade, considerable 
work has been done with the most realistic nuclear potentials to determine the magnitude 



of the energy gap in the triplet superfluid neutron matter for different densities jsj- 10|. So- 
phisticated calculations have shown that, besides the above one-component state, there are 
also multicomponent states involving several magnetic quantum numbers that compete in 
energy and represent various phase states of the condensate dependent on the temperature. 

Whether the phase transitions modify the spectrum of low-energy excitations and the 
intensity of neutrino emission from the volume of neutron stars is the question we try to 
answer in this paper. Theoretical investigation of low-energy excitations responsible for the 
neutrino emission by the neutron triplet superfluid liquid is conducted first. Until recently, 
the only known excitations able to decay into neutrino pairs were the broken pairs. It is well 
known that the neutrino emission caused by the pair-recombination processes in the neutron 



triplet superfluid liquid can dominate in the long-term cooling of neutron stars Uj . We will 
consider also the collective excitations in the timelike domain of energies and momenta, 
which can also be responsible for the intense neutrino emission. Since the neutrino emission 



in the vector channel of weak interactions is strongly suppressed [12| we will focus on the 
collective spin-density oscillations that can decay into neutrino pairs through neutral weak 
currents. 
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Previously spin modes have been studied in the p-wave superfluid liquid ^He 
The pairing interaction in ^He is invariant with respect to rotation of spin and orbital 
coordinates separately. In this case, the spin fluctuations are independent of the orbital 
coordinates. In contrast, the triplet-spin neutron condensate arises in high-density neutron 
matter owing mostly to spin-orbit interactions that do not possess the above symmetry. 
Therefore the results obtained for liquid ^He cannot be applied directly to the superfluid 
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neutron liquid. 

Recently spin waves with the excitation energy smaller than the su per fluid energy gap 
were predicted to exist in the superfluid condensate of neutrons 12[. The neutrino 
decay of such spin waves is important for thermal evolution of neutron stars with 
the conventional one-component ground state with rrij = 0. In this paper, we consider 
spin-density excitations for the other superfluid phases, which can be preferred at some 
temperatures. 

We will not consider the spin oscillations of the normal component. These soundlike 
waves that transfer into the ordinary spin waves in the normal Fermi liquid above the critical 
temperature cannot kinematically decay into neutrino pairs. Instead, we will focus on the 
spin excitations of the order parameter, which are separated by some energy interval from 
the ground state and are kinematically able to decay into neutrino pairs. The dispersion 
equation for such waves in the superfluid one-component condensate with mj = was 



derived in Ref. 



12| in the BCS approximation. In this paper we study the collective 



spin excitations in multicomponent phases of the condensate, while taking into account the 
Fermi-liquid interactions. 

This paper is organized as follows. Section II contains some preliminary notes and outlines 
some of the important properties of the Green functions and the one-loop integrals used 
below. In Sec. Ill we discuss the renormalization procedure which transforms the standard 
gap equation to a very simple form valid near the Fermi surface. In Sec. IV we derive the 
effective ordinary and anomalous three-point vertices responsible for the interaction of the 
multicomponent neutron superfluid liquid with an external axial-vector field. We analyze 
the poles of anomalous vertices in order to derive the dispersion of spin-density oscillations 
in the condensate. In Sec. V we derive the linear response of the multicomponent superfluid 
neutron liquid onto an external axial-vector fleld. In Sec. VI we briefly discuss the general 
expression that relates the neutrino energy losses through neutral weak currents to the 
imaginary part of response functions. We derive the neutrino losses caused by recombination 
of broken Cooper pairs and by decay of spin waves. Finally, in Sec. VII, we evaluate 
neutrino energy losses in the multicomponent superfluid neutron liquid undergoing the phase 
transition. Section VIII contains a short summary of our flndings and the conclusion. 

Throughout this paper, we use the standard model of weak interactions, the system of 
units h = c = 1 and the Boltzmann constant /c^ = 1- 
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II. PRELIMINARY NOTES AND NOTATION 



The spin-orbit interaction between quasiparticles is known to dominate in the nucleon 
matter of high density. The most attractive channel corresponds to spin, orbital, and total 
angular momenta s = 1, / = 1, and j = 2, respectively, and pairs quasiparticles into the 
states with rrij = 0, ±1, ±2. The substantially smaller tensor interactions lift the strong 
paramagnetic degeneracy inherent in pure pairing and mix states of different magnetic 



quantum numbers Sj- 10|]. The admixture of the '^^2 state, which arises because of the tensor 



interactions, is known to be small and does not affect noticeably the excitation spectra [19|. 
Accordingly, throughout this paper, we neglect small tensor forces but consider the case 
of pairing into the multicomponent (mj-mixed) states corresponding to the phases of the 
realistic superfluid condensate. The pairing interaction, in the most attractive channel, can 
then be written as [l| 

Q^ap,^s (P, P') = V {p, p') ^ [hm^ {ri)crg) [gdrh^^^ '^'^'O ^5 ' 

where V {p,p') is the corresponding interaction amplitude, g = ppM* /n'^ is the density of 
states near the Fermi surface, <t = (ai, d'2, 0-3) are Pauli spin matrices, g = id'2, and b^^. (n) 
are vectors in spin space that generate the standard spin-angle matrices, so that 

hm^{n)a-g = ^ Q^a/3|lm J (llmsm;|2mj) Yi,^, (n) . (2) 

ms+mi=mj ^ ' 

These are given by 



bo = VT72(-ni, -n2,2n3) , bi = - v^S/i (ns, mg, ni + mz) , 

ba = V374 {ni + in2, mi - n2, 0) , b_m, = h*^^ , (3) 

where ni = sin ^ cosy), n2 = sin^siny^, and na = cos^. The vectors are mutually orthogonal 
and are normalized by the condition 

^h*^>hm, = Sm,m'. (4) 

The triplet order parameter D = Dap (n) in the neutron superfluid represents a symmetric 
matrix in spin space =t, i), which can be written as 

D (n) = ^ A^^, (6-b„J ^ . (5) 



We are mostly interested in the values of quasiparticle momenta p near the Fermi surface, 
p ~ Pp^ where the partial gap amplitudes A^^. {p) Amj (pf) are almost constants, and 
the angular dependence of the order parameter is represented by the unit vector n = p/p, 
which defines the polar angles {6, if) on the Fermi surface. 

The ground state (|5]) occurring in neutron matter has a relatively simple structure (uni- 
tary triplet) 

(n) = A b (n) , (6) 

where A is a complex constant (on the Fermi surface), and b (n) is a real vector which we 
normalize by the condition 

r rlr\ 

1. (7) 



Various sets of the gap amplitudes A^^ in Eq. correspond to the various phases of the 
condensate considered further. 

By making use of the adopted graphical notation for the ordinary and anomalous prop- 
agators, G = — » — , G^{p) = — • — , -F*-^-* = * * , and F^^-' = « » , we employ the 
Matsubara calculation technique. Then the analytic form of the propagators is as follows 



20|, |21| 



G {rin,p) = G {r]n,p)6ap , G {r]n,p) = G (?7„,p)(5«^, 

(r^„, p) = F (t]^, p) h&g , F^^) (r^„, p) = F (r/„, p) g&h , (8) 

where the scalar Green functions are of the form G~ (?/„, p) = G {—r]n, — p) and 

Vn + EI vl + El 

Here, ?7„ = m {2n + 1)T with n = 0,±1,±2... is the fermionic Matsubara frequency and 
= vf {p — Pf) with vf the Fermi velocity. The quasiparticle energy is given by 

El=el + (n) , (10) 

where the (temperature-dependent) energy gap Ab (n) is anisotropic. In the absence of 
external fields, the gap amplitude A (T) is real. 

In general, the Green functions ([8]) should involve the renormalization factor a ~ 1 



independent of a;,q, T (see e.g., |2l|). The final outcomes are independent of this factor; 



therefore, to shorten the equations, we will drop the renormalization factor by assuming 
that all the necessary physical values are properly renormalized. 

Finally we introduce the following notation used below. We designate Xxx' ij^i n, q; T) as 
the analytical continuations onto the upper-half plane of complex variable u of the following 
Matsubara sums: 

(w„^, n, q; T) = T ^ ^ y c?£:pX (^r^^ + w^, p+| j X' (^?7„, p-|j . (11) 

where X, X' G G, -F, G~ and Um = ^inTm with m = 0, ±1, ±2, ....These are functions of w, 
q, and the direction of the quasiparticle momentum p = pn. 

The loop integrals fill I) possess the following properties, which can be verified by a 
straightforward calculation (the same relations have been obtained in Ref. {22! for the 
case of singlet-spin condensation): 

^G-G = '^GG- ! ^GF = —^FG ; '^G- F = —^FG- ; (12) 

^G-F + ^FG = ' (1^) 

Tq-p — IpG = —^IpF ■ (14) 

For arbitrary w, q, T one can also obtain 

Igg- + b'iFF = A+ "^^l^^Ipp , (15) 

where v =VF'n, and 

^(n) ^ [Xa-G(n) +6Mn)X^p(n)]_^^^^^ . (16) 

In the case of a triplet superfluid, the key role in the response theory belongs to the 
loop integrals Xpp and [Igg i b'^^pp) . For further usage we indicate the properties of these 
functions in the case of u > and q — i- 0. A straightforward calculation yields 

hmipp = I {n,uj) = — — -2tanh— , (17) 

and 

{Igg + PIff)^^, = 0, (18) 

{iGG-b'ipp) =-2Pl. (19) 
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III. GAP EQUATION 

The standard gap equation [3| involve integration over the regions far from the Fermi 
surface. This integration can be ehminated by means of the renormahzation of the pairing 
interaction 22|]. We define 



T) = V {p,p') - g-' I {p,p") (GG")'; {p",p'; T) , 



(20) 



where the loop {GG )^ is evaluated in the normal (nonsuperfluid) state. Then it can be 



shown 12| that we may everywhere substitute V^"^^ for V provided that at the same time, 
we understand by the GG~ element, the subtracted quantity GG^ — {GG^)j^ [{GG~)j^ is 
to be evaluated for a; = 0, q = in all cases]. 



The function (fT6|) is now to be understood as 

A (n) ^ [la-a - T^c-g)^ + b'lj.^]^^,^^^, (21) 

and the standard gap equations can be reduced to the form 

^m, = -Al^M I ^h*^^{n)h{n)A (n) , (22) 

which is valid in the narrow vicinity of the Fermi surface where the smooth functions A.^^ (p), 
{p,p'), and A(p) may be replaced with constants A (p) ^ A (pp) = A, etc. 
The function fl2ip can be found explicitly after performing the Matsubara summation: 

... r . f 1 , + A262 1 e\ 

A(n) = - d£ { , tanh tanh — . (23) 

^ ' 2 70 \^e^ + /\W 2T e 2T) ^ ^ 

IV. EFFECTIVE VERTICES 

The field interaction with a superfluid liquid should be described with the aid of four 
effective three-point vertices. There are two ordinary vertices, T(n) , -f^ (n) = t^(— n), 
corresponding to creation of a particle and a hole by the field (which differ by the direction of 
fermion lines), and two anomalous vertices, T^^^ (n) and T*^^^ (n), corresponding to creation 
of two particles or two holes. 

The anomalous effective vertices are given by infinite sums of the diagrams, taking into 



account the pairing interaction in the ladder approximation 23|. The ordinary effective 



vertices incorporating the particle-hole interactions can be evaluated in the random-phase 
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approximation |2J]. This can be expressed by the set of Dyson equations symbohcally 
depicted by graphs in Fig. [1] 





>■ = >■■ ^ x> - ^ x> - ' O - ' x>- 

^^^^ ^CZ-^ ^^z^ ^xz^ 



FIG. 1. Dyson equations for full ordinary and anomalous vertices. The particle-hole interaction is 
shown by the shaded rectangle. Wavy lines represent the pairing interaction. 

In these diagrams the shaded circle is the full ordinary vertex, and the shaded triangle 
represents the anomalous vertex. The particle-hole interaction is shown by the shaded 
rectangle. Wavy lines represent the pairing interaction. The first diagram on the right-hand 
side of the first line is the three-point vertex of a free particle. 

In our analysis, we shall use the fact that the Fermi-liquid interactions do not interfere 
with the pairing phenomenon if approximate hole-particle symmetry is maintained in the 
system, i.e., the Fermi- liquid interactions remain unchanged upon pairing. Since we are 
interested in values of quasiparticle momenta near the Fermi surface, p ~ p^n, the Fermi- 
liquid effects are reduced to the standard particle-hole interactions: 

Qda^,l35 (nn) = f (nn ) 6af)S^5 + g (nn) CTapcr-yS ■ 

We are interested in excitations able to decay into neutrino pairs through neutral weak 
currents. Since the neutrino emission in the vector channel of weak interactions is strongly 
suppressed in nonrelativistic media {2^, 2^], [l2|, we will focus on the interaction of the 
superfiuid Fermi liquid with an external axial-vector field. In the nonrelativistic case, the 
bare axial-vector vertex is given by the spin matrices &. (We neglect a small temporal 
component that arises as the relativistic correction.) 

After the proper renormalization of the pairing interaction the equations for the axial- 
vector vertices can be reduced to the following analytic form (for brevity, we omit the 



dependence of functions on u and q ): 

d-a! 

T(n) = & + & 



-S(nn') XgfTi (tT (6-b) 



(ab) 



-XccTr [erf] + XffTi [ct (o-b) gf g {crh)] , 



(24) 



t«(n) = J]6-b„^(n)^\/M / ^^Xc,G-Tr 
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X^^Tr[(ab:;^.) (6-b) ^t(2) (6-b) 



XgfTi 



(o-b;;^) T (ab)] +x^G-Tr [(o-b;;,^,) (^b) ^r-^] , (25) 



rfn' 



t(2) (n) = '£gahl^^{n)V^^) / {x,,-^Tr [(^b„J ^t^^) 



-X^^Tr [((Tb^J (6-b) (6-b) 
+XG-FTr [(o-b^J ^f"^ (6-b)] -X^^cTr [(<Tb„J (ctB) t] 

Inspection of the equations reveals that the solution should be of the form 

f(n) = (n) 6- , gf- (n) g = (p (-n) a 

t« (n,c.) = 5^BW(c.) {ah^^)g, 

t(2) (n,c.) = ^B(^)(a;)^(^b: 



(26) 



(27) 
(28) 

(29) 



After this substitution and summation over spins, one can obtain a set of equations for 



i±(n) = -(0(n)±0(-n)) 



(30) 



and 



B^, = ^(b«±(-)"^bL1^,) . (31) 
The application of a little algebra using Eqs. (IT^ -( IT^ results in the following equations: 

0+ (n) = 1 + 1 ^0 (nnO {Igg - h'Tpp)^, 0+ (n') 

^0(nn')Xp^(n')^(b^^.xb)^, , (32) 



(n) = / ^0 (nn') (Xgc + 0- (n') 



(33) 



m'- 



2(b:;,^xbjx.. 



X 



(34) 



dn 



— (qv)^ 

2A2 



2P b* >^. +2 b* ,.b bb 



X 







(35) 



In obtaining the last two equations, we used the gap equation fl22|) and the identity f[T5|) . 

Further simphfications are possible due to the fact that in Eq. (135|) do not couple to 
external fields. Even if the eigenoscillations of Bl, exist it is unclear how this mode could 
be excited. Therefore one may assume that Eq. ( !35ll has only the trivial solution = 0. 
This simplifies Eq. ( l32l) which is now uncoupled. 

The amplitudes of Fermi-liquid interactions can be expanded into Legendre polynomials 
and written in terms of an infinite set of Landau parameters. In the axial channel, this gives 



0(nn') = X^0;P;(nn') . 



(36) 



1=0 



We now expand the functions (j)± (n) over spherical harmonics (n). It is apparent that 
the function (n) = 0+ (— n) contains only even harmonics, 

I 

0+ (n; U), q) = 0+„ (w, q) yj,„ (n) , (37) 

i=even m=—l 

while 0_ (n) = — 0_ (— n) consists of odd harmonics, 

(38) 



(n; u, q) = ^ ^ 0, „ (w, q) F/,^ (n) 



l=odd m=—l 



An 



Making use of the relation 

dnYl^ (n) g (nn') = ^S'^a,;!^,; (n') 

which follows from the expansion fl36|) . we arrive at the final set of of equations, 

/'=even m' 



Lm 



21 + 



Y E E 'i'v.n' j dnY*m i^GG + b'lpp) Y, 

l'=odd ml 



^ /ATTgi . qvF 
+ 2/ + l' A 



Y^n.J^ (cOS^qn) (b^^. Xb) IffY* 



(39) 



(40) 



m'. 



dn 



^mj ,m 



A 



bl b A 



+ 



An 



2A2 



A 



-b:„^,b^-2fb:,^b) (bb^. ) IX, 



PF 



;=even m 

7 ™ 



(41) 

=odd m 

where is the angle between the transferred momentum and the direction of quasiparticle 
motion. 

Since / can take all values from zero to infinity, a general solution cannot be given in 

closed form. As in the case of a normal Fermi-liquid, a closed solution may be obtained if 

we set = for / > 1. We adopt this approximation and consider the solutions with q = 0. 

With the aid of Eqs. (HZD-CH]), we find 0,;^ = and 

1 



1 + 2go (6^X(n; w)) 

Hereafter the angle brackets denote angle averages, (...) = {An)~^ j dn.... Inserting these 
functions in Eq. (HTl) we obtain the set of equations for B+ , 



bl ,b^. U) -5. 



A 



•'A 



2A2 



blb^/ 1 X 



2( (bl b) (bb^, )X 



^^^(fc xb)X 

1 + 2go (&2X) A \ V ™^ 



(43) 
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The unit vector b used here was defined by Eq. ([6]). 

The exphcit evaluation of Eq. ( H3l) for arbitrary values of u and T appears to require 
numerical computation. However, we can get a clear idea of the behavior of this function 
using the angle- averaged energy gap A^6^ — > (A^6^) = in the quasiparticle energy (fTOj) . 
(Replacing the angle-dependent gap in the quasiparticle energy by its average has been 
found to be a good approximation ^, 12|, 19|.) In this approximation, the functions 
X {uj, T; n) — Xav (w, T) and A (T), in Eqs. f H2|) and ( l43l) . can be moved beyond the angle 
integrals. With the aid of normalization condition (^^, we find 

1 



0ao (w> T) 



Using also the fact that 



1 * 1 



5. 



rrijm- > 



A 



(44) 



(45) 



and substituting 



1 



from Eq. 



1 + 2goXa, 

3]), we obtain the set of linear equations for Bm^. with mj = 0, ±1, ±2, 

2 



4A2 



Br 



bbr 



(46) 



For further progress, we need to define a particular form of vector b that characterizes the 
ground state of the condensate. The general form of a unitary state is to be written as 



Ao 
A 



A 



(bi - b_i: 



A2 

A 



(b2 + b_2) 



(47) 



with 



A^ = A^ + 2Al + 2Al . (48) 

By utilizing notation adopted in Refs. js], jo], where Ai = -\/6Ai/Ao and A2 = -\/6A2/Ao, 
from Eq. (H7I) we obtain the general form of the properly normalized vector b: 



lAo 
2 A 



--r- -Til + niX2 - n-iXi , -n2 - n2\2 , 2^3 



The solution to the set of linear equations ( H6|) is found to be 

Ai 



Bn 



1 

2 A3cj2 
11 



-I (0,^,0) 



(49) 



(50) 



where u = u/ (2A), and we use the following notation: 



i = ^§ (5 + 2X1 - 2A2 + + ^(1 + A^)^ (4Af + (A2 - 3)^) 
I (^5 + 2X\ - 2X2 + Xl- ^(1 + A2)' (4A? + (A2 - 3)') 



40 A2 



(53) 
(54) 
(55) 



,_Al fl Xi2 + X2 



A2 \5 10 3+ A2 
20 A 



^4 = 7^^(A? + 2A2 + 2A2 + 4) 



^^(il-X2? + Sil + SX2 



20 A2 



2A, 



By taking into account that B^. = 0, we have B 

1 



(1) 



(2) 



and 



t(2) fn, 



1 + 2(7oXa. 

1 



CO 



^Bm^ {uj) (o-b^J ^ , 

rrij 



(56) 
(57) 



1 + 2(?oXa, 

In the last equation, the property ([3]) is used. 

The full ordinary axial-vector vertices can also be simplified because = and, there- 
fore, 0(— n) = 0(n). From Eqs. (!37|) . ( H2|) . and (14^ . we find (n; w, q) = 0^Q(cj,q) and, 
thus. 



1 



1 



1 + 



-cr , r n,cj 



According to modern theories 



1 + 20oXa, 



10l |. there are several multicomponent states that com- 



pete in energy depending on the temperature. Accordingly the phase transitions occur 
between these states when the temperature goes down. The possible phase states of the 
F2 condensate are cataloged in Ref. {q]. Immediately below the critical temperature, 
the superfluid condensate can appear either in one of the two two-component phases, 

Ao 1 



0±3 : ^ = - , Ai = , A2 = ±3 



(58) 
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or in the one-component phase, 

mj = 0: Ao = 1, Ai = 0, As = . (59) 

These lowest-energy states are nearly degenerate. The higher group is composed of the 
phases 



Oi-- x= ^ / ' A, = ' /2(l7-3V2l), A. = |(v/21-4) (60) 



and 

02 : 4^ = _ / ^ , Xi = ^J2 (l7 + 3v^), As = (v^ + 4) . (61) 



According to the above calculations the energy split between the two groups shrinks along 
with the temperature decrease and results in the phase transition at T = 0.7Tc. 

The effective vertices and the polarization tensors for each of the above phases can be 
obtained with the aid of Eqs. (!50|) -( [52|) and (!67|) . We found that for all of the above- 
mentioned phases, the anomalous vertices can be written universally in the form 

''- = (x)';5^K''™,-s) (62) 



where 



(b*xb) = --i=^(o,v^Ai,0 



2V6 A 

(hl.xh) = H (3 + As) , T (3 - As) , -^X^) , 

(bisxb) = (T^Ai, -Ai, t2^As) . (63) 

Poles of the anomalous vertex indicate eigenmodes of the order parameter. The pole at 
o)^ = 1/20 signals the existence of undamped collective spin oscillations of the energy 

a;, (g = 0) = A . (64) 

Three important conclusions follow immediately from this simple formula: (a) The spin- 
density waves are of the identical excitation energy in all phases of the superfluid neutron 
liquid, (b) Fermi-liquid interactions do not influence the spin-wave energy in the condensate, 
(c) The spin waves are kinematically able to decay into neutrino pairs through neutral weak 
currents. 
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V. POLARIZATION FUNCTIONS. 



Unfortunately, the Landau parameter go for the particle-hole interactions in asymmetric 
nuclear matter is unknown. Therefore, in evaluating neutrino energy losses, we simply 
neglect the Fermi-liquid effects by taking go — and consider the axial polarization in the 
BCS approximation. The latter can be obtained in the form 12| 



2A 



XppTr 



diff^ (o-b) - CTi (o-b) gf^ 



(2) 



(65) 



where z, j = 1, 2, 3 and the anomalous axial-vector vertices T*^^'^^ are given by Eqs. f l56|) and 
( l57l) . As in the above, we focus on the case q = and omit for brevity the dependence on 
n and u. By using Eq. f ITI?]) and applying a little algebra, we obtain 



-Aq{{5%^ -W)X) 



+ 



+ 20) - AV5 A2 



The pole location on the complex u plane is chosen to obtain the retarded polarization 
tensor. 

For further evaluation of the polarization function we will again use the angle-average 
approximation by replacing the angle-dependent gap in the quasiparticle energy ( TTUj) with 
its angle average, 6^A^ -> (6^A^) = A^. Then the function 



2^av (c^) 



de 



A^ 



^ tanh — 



(66) 



/o E 4^2 _ + ^0)^ 

with E"^ = e"^ + 4A^, is isotropic and can be moved beyond the angle integral. In this 
approximation, we obtain 



+ 



A 



;cc; + zO)' - AV5 '"'"'A 



(67) 



Below we use the retarded polarization tensor for calculation of the neutrino energy losses 
from superfiuid bulk matter of neutron stars. In this calculation one can neglect the temporal 



and mixed components of the tensor occurring as small relativistic corrections 



25| 



14 



VI. NEUTRINO ENERGY LOSSES 



We will examine the neutrino energy losses in the standard model of weak interactions. 
Then after integration over the phase volume of freely escaping neutrinos and antineutrinos 
the total energy which is emitted per unit volume and time can be obtained in the form (see 



details, e.g., in Ref. |28|) 

' " 192^^ X J ^ exp(^) -1 ^ ^ " ~ ^^"^ ' ^ ^ 

where Gp is the Fermi coupling constant, Ca = 1-26 is the axial-vector weak coupling 
constant of neutrons, AC = 3 is the number of neutrino flavors, O (x) is the Heaviside step 
function, and fc^ = (w, q) is the total energy and momentum of the freely escaping neutrino 
pair (/i, = 0, 1, 2, 3). 

In Eq. ( I68ll . we have neglected the neutrino emission in the vector channel, which is 
strongly suppressed due to conservation of the vector current. Therefore the energy losses 
are connected to the imaginary part of the retarded polarization tensor in the axial chan- 
nel, n^'^ ~ 6'^^6''^ ImlVp^. The latter is caused by the pair breaking and formation (PBF) 
processes and by the spin- wave decays (SWDs). These processes operate in different kine- 
matical domains, so that the imaginary part of the polarization tensor consists of two clearly 
distinguishable contributions, Imll^ = Imllpgp + Imllg^j-), which we will now consider. 

A. PBF channel 

The imaginary part of Xav, which arises from the poles of the integrand in Eq. ( 166|) at 
uj = ±2E is given by 

ImXav (UJ>0) = , ^ tanh . 

The PBF processes occur if a; > 2 A, while the energy fl64l) of the spin waves is considerably 
smaller. Therefore we may neglect the term A^/S relative to uP' in the denominator of Eq. 
(EH). Then, 



X ; ^ tanh — , (69) 

uVoo^ - 4A2 4T ' ^ ^ 
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where 



The tensors {bibj") and can be found with the aid of Eqs. ( H9|) and (l63l) . A straightforward 
calculation gives 



[bibi) 



1A| 
6 A2 



1 + 2A2 + 

-A1-A1A2 4 + XI J 



(70) 



and 



1 

A.. = _flo 
12 A2 



A? + (A2 + 3)' 







3A 



l(A2 + l)\ 



V 



(71) 



4A2 + (A2 - 3)' 

3Ai(A2 + l) XI + 4X1 J 

Inserting the imaginary part of the polarization tensor into Eq. ( |68|) . we calculate the 

contraction of Imllpgp with the symmetric tensor k^^k,^ — k'^g^y. This gives 

1 . , ..... r dco eitu^- 4A2) 

'0 



967r' 



{l + exp^Y Vco^ - 4A2 
r f - - 3 A^ A^ 



(72) 



Integration over d^q can be done in cylindrical frame, where qi = g^cos^, q2 = g_Lsin$, 
and qs = qz- This results in the neutrino energy losses in the form 

u^du A2 



eOvr^ V" 4AV -------- J^^ (i + expi^)Vu;2-4A2 

In obtaining this expression, the following fact is used: 



(73) 



Tr Aij = Ti (bibi) 



A2 n 



xl + lxl + i 



1 . 



3 A2 \3 ' ' 3 

With the aid of the change u = 2Ta/x^ + A^/T^, one can recast Eq. (!73|) into the form: 

2 f . „ An \ . - , of ^ dx 



CPBF 



ISvrS 



4-3||) GIcIKpfM*TV f - 



(1 + exp z) 



(74) 



where z = + and y = A (T) /T. 

For a practical usage from Eq. fTMl) . we find 



6PBF = 5. 85 X 10^° (^) T.W.ClFpBF 



(y) 



erg 
crrfis 



(75) 
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where M and M* = pp/vp are the bare and effective nucleon masses, respectively; Tg 
T/IO^K, and 



^4 

^1 + exp z 

The neutrino energy losses, as given by Eq. (!74|) with good accuracy reproduce the result 



PBF (y) = ( 4 - 3^) rfx^^^l-^. (76) 



obtained in Ref. 



12] for the one-component phase rnj = 0, where Aq/A = 1. It is necessary 



to notice that Eq. (!76|) obtained in the angle-average approximation is much simpler for 
numerical evaluation than the "exact" expression which contains additionally the angle 
integration [l^. To avoid possible misunderstanding we stress that the gap amplitude 
A (T) in Eq. fITBl) is \/2 times larger than the gap amplitude Aykl used in Ref. ll| , where 



the same anisotropic gap An = A6 (n) is written in the form An = AyklV^ + 3 cos^ 9 = 
AYKiV^b^n). In other words, (A^) = A^ = 2A^^^. 

The small difference in the gap amplitudes, 5A^/A^ ~ 2%, inherent for various phases 
of the condensate is crucial for the phase transitions [lOj, but this small inequality can be 
disregarded in evaluation of the neutrino energy losses. Therefore the efficiency of PBF 
processes in various phases of the superfluid condensate is proportional to the factor 

3 A2 

ft:pBF = 1-4^ 

For the one-component phase rrij = one has kpbf = 1/4; in the case of 0±3 condensates, 
KPBF = 13/16. For the d phase, we obtain kpbf = (l73 - 9721) /224 ~ 0.5882, and for 
the O2 phase, we have kpbf = (l73 + 9^21) /224 ~ 0.95644. 



B. SWD channel 

In the frequency domain < 00 < Ab, the imaginary part of the weak polarization tensor 
( |67|) arises from the pole of the denominator at a; = A/-\/5 and can be written as 

ImnSWD (c^ > 0) = -^nV5Ag5 (a; - A/Vtj I,, (uj) ^A,, . (77) 

Inserting this expression into Eq. (1^5]) and performing trivial calculations, we find 

GlClK ,,,Al {A/V^y f^de A2 ^ K 
= At ,,p(AMj_ J o E E^-aW ^''2T- ^''^ 
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One can neglect in comparison to E"^ > 4A^ in the denominator of the integrand. 

With this simphfication, we obtain 



7 



-wo - i^gVMT' ^^^^^;^_^ (^) X. , (79) 

where y = A (T) /T, and 

Xo (y) = / — — tanh . (80) 

This expression can be recast into the traditional form 

.swo = 2. 74 X 10" ( ^) (^) rXCi§M^ ^ . (81) 

\M J \McJ exp (y/v5) - 1 cm^-s 

From this expression, it is seen that the efficieency of spin-wave decays is proportional to 
i^swD = Al/A"^. For the one-component condensate rrij = 0, kswd = 1; in the case of 0±3 
condensates, kswd = 1/4. For the Oi phase, we obtain kswd = (l7 + 3a/21) /56 ~ 0.549 07, 
and for the O2 phase, /tgwD = (l7 - 3^21) /56 ~ 5.81 x lO^^. 

VII. PHASE TRANSITIONS AND EFFICIENCY OF NEUTRINO EMISSION 

For numerical evaluation of the neutrino losses, it is necessary to know the function 
y = A (T) /T, which in general is to be found with the aid of gap equations. However, as 
mentioned above, the difference in the gap amplitudes for various phases can be neglected in 
evaluation of the neutrino energy losses. This substantially simplifies the problem because for 
the case mj = the function is well investigated. We can adjust, for example, the simple fit to 



Aykl (T) it = Mb (t), as suggested in Ref. [ll|, where r = T/T^. Taking into account that, 
in Ref. [Uj, the gap amplitude Aykl (T) is defined by the relation = Ayj^^ (1 + 3 cos^ 6), 
while our definition is A^ = ^A^ (1 + 3 cos^ 9), we obtain y (r) = \/2\ib (t)- 

In Fig. [2] we compare the PBF and SWD neutrino emissivity for various phases of 
superfiuid neutron matter. The temperature dependence of the emissivity is evaluated at 
Pp = 2.1/m~^. We set the effective nucleon masses M* = 0.7M; the critical temperature 
for neutron pairing is chosen to be = 3 x lO^i^'. 

One can see that the decay of spin waves into neutrino pairs is very effective at low 
temperatures, when other known mechanisms of neutrino energy losses in the bulk neutron 
matter are strongly suppressed by superfluidity. Maximal neutrino emission in the SWD 
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FIG. 2. Temperature dependence of the neutrino emissivity due to recombination of Cooper pairs 
(PBF) and due to decay of spin waves (SWD) for pp = 2.1 fm^^ and Tc = 3 x 10^ K. The effective 
mass is taken to be M* = 0.7M. The curves of different styles correspond to various phases of 
superfluid condensate, which are discussed in the text. 

channel occurs in the one- component phase. In contrast, efficiency of the PBF channel is 
maximal in the O2 phase. 

In Fig. [3] we demonstrate the total neutrino emissivity e = epBP + cswd versus the 
temperature by assuming that the phase transition occurs at T = 0.7Tc. The phase transition 
(if it occurs) leads to a sharp increase in the neutrino energy losses followed by a decrease, 
along with a decrease in the temperature that takes place more rapidly than it would without 
the phase transition. 

According to the minimal cooling paradigm [gJ, [3], along with lowering of the temper- 
ature, the star continues to lose its energy by radiating low-energy neutrinos via the PBF 
processes untill a photon-cooling epoch enters at T ~ O.lTc. At this latest stage of the 
cooling, all mechanisms of neutrino emission from the inner core are suppressed greatly by 
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FIG. 3. Total energy losses e = epBF + eswD versus the temperature. The superfluid condensate 
undergoes the phase transition at T = 0.7Tc. Parameters of the medium are the same as in Fig. 2. 

the neutron and proton superfluidity, and the 7 radiation from the star surface is considered 
as the main mechanism of the star coohng. 

Given the strong dependence of the PBF and SWD neutrino emission on the temperature 
and the density, the overall effect of the SWD processes can only be assessed by complete 
calculations of the neutron star cooling which are beyond the scope of this paper. A rough 
estimate can be made by considering a simplified model of the superfluid core of the density 
2po enclosed in the volume 7 x 10^^ cm? . In Fig. IHwe demonstrate the total neutrino energy 
losses caused by PBF and SWD neutrino emission from the star volume in comparison with 
the surface photon radiation. The latter is taken as in Fig. 20 of Ref. j^l. There is a 
relatively large range of predicted values for T^; therefore, we show the PBF and SWD 
neutrino luminosities for Tc = 3 x 10^ K and for Tc = 3 x 10^ K. This simple estimate shows 
that the neutrino emission caused by spin-wave decay (SWD) can dominate the 7 radiation 
within some temperature range, which was previously considered as the photon-cooling era. 
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FIG. 4. Temperature dependence of the neutrino luminosities from the model neutron star due to 
recombination of Cooper pairs (PBF) and due to decay of spin waves (SWD) for pp = 2.1 fm"^. 
The effective mass is taken to be M* = 0.7A/. Solid lines correspond to Tc = 3 x 10^ K; dot-dashed 
lines correspond to Tc = 3 x 10^ K. Volume of the triplet condensate is estimated as 7 x 10^^ cm^ . 
The dashed line is the energy losses per unit of time due to surface 7 radiation, as calculated in 
Ref. y. 



To get an idea of the lowering of the cooling trajectory due to SWD neutrino emission 
we consider a simple model of cooling of the above core of the density 2po and of the volume 
7 X 10^^ cm^ enclosed in a thin envelope typical for real neutron stars by assuming that the 
surface temperature is connected to the central temperature by the formula in Ref. 29| . We 
assume also that the bulk matter consists mostly of superfluid neutrons with = 
and contains a small admixture of normal (nonsuperfluid) protons and electrons (the proton 
fraction Xp = 0.05), so that the total specific heat Cy consists of the three corresponding 
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FIG. 5. Transition of cooling trajectories between a model with a surface photon radiation and 
PBF neutrino emission and a model including a surface photon radiation and PBF+SWD neutrino 
processes. Tc = 3 x 10^ K. Details of the calculation are described in the text. 



contributions, as described in Ref. (30|- Under these conditions, the coohng equation 

-Q , (82) 



dT 



can be solved numerically with T = at the initial moment. We have solved this equation 
for the case when Q = + Qpbf and for the case when Q = + Qpbf + Qswd- The 
result is shown in Fig. [5l 



VIII. SUMMARY AND CONCLUSION 



Let us summarize our results. We have studied the linear response of the superfluid 
neutron liquid to an external axial-vector field. The calculation is made for the case of 
a multicomponent condensate involving several magnetic quantum numbers and allows us 
to consider various phases of superfluid neutron liquid. In order to estimate the neutrino 
energy losses, while taking into account possible phase transitions, we have considered the 
low-energy excitations of the multicomponent condensate. 

Along with the well-known excitations in the form of broken Cooper pairs, we consider the 
collective waves of spin density, which are known to exist in the one-component condensate 
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at very low energy [18|] . Our theoretical analysis predicts the existence of such waves in all of 
the multicomponent phases we have considered. We found that the excitation energy of spin 
waves is identical for all of the phases and is independent of the Fermi-liquid interactions. 
In the angle-average approximation, the energy of spin-density oscillations is estimated as 
(g = 0) ~ A/y/5. 

Neutrino energy losses caused by the pair recombination and spin-wave decays are given 
by Eqs. (!7il) and ( 1791) . respectively. Because of a rather small excitation energy, the decay 
of spin waves leads to a substantial neutrino emission at the lowest temperatures T ^ T^, 
when all other mechanisms of the neutrino energy losses are killed by a superfluidity. We 
have evaluated the neutrino energy losses for all of the multicomponent phases that might 
represent the ground state of the condensate according to modern theories. 

Finally we have evaluated the temperature dependence of neutrino energy losses from the 
superfiuid neutron liquid in the case when the phase transition occurs in the condensate at 
the temperature T = O.TTc. Our estimate predicts a sharp increase of the neutrino energy 
losses followed by a decrease, along with a decrease of the temperature that takes place more 
rapidly than it would without the phase transition. 

Since the neutron triplet-spin pairing occurs in the core which contains more than 90% of 
the neutron star volume, the neutrino processes discussed here could influence the evolution 
of neutron stars. 
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